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ABSTRACT 



This work considers flows from an accretion disk corotating with the aligned dipole magnetic field of a 
rotating star. Ideal magnetohydrodynamics (MHD) is assumed with the pressure and density related as 
p cx p 1 and with pv 2 <C B 2 /47r, where v is the flow velocity. This limit corresponds to the Alfven radius 
for the disk accretion larger than the corotation radius and is of particular interest because flow can be 
solved rigorously. Transonic flows, which go from subsonic motion near the disk to supersonic inflow near 
the star, are shown to be possible only for a narrow range of R d ~ r c = (GM/fi 2 ) 1 / 3 , where R d is the 
radius at which the dipole field line intersects the disk, r c the corotation radius, M the mass of the star, 
and CI its angular rotation rate. Over a larger range of Rd < r c , subsonic flows from the disk to the star 
are possible. The transonic flows have very different behaviors for 7 > 7/5 and < 7/5. In both cases, 
the plasma flow velocity v (which is parallel to B) increases with decreasing distance R from the star. 
However, for 7 > 7/5, the Mach number M. = |v|/c s (with c s the sound speed) initially increases with 
R decreasing from Rd, but for R decreasing from w 0.22Rd (for 7 = 5/3) the Mach number surprisingly 
decreases. In the other limit, 7 < 7/5, A4 increases monotonically with decreasing R. Application of 
these results is made to funnel flows to rotating magnetized neutron stars and young stellar objects. The 
spatial variation of both the flow velocity and the Mach number are crucial to the nature of the standing 
shock wave near the star and to the determination of emission line profiles. 

Subject headings: physical processes: accretion, accretion disks: hydrodynamics: plasmas — stars: 
formation: magnetic fields: pulsars 



an aligned rotator, but these authors do not analyze the 
Bernoulli equation and do not consider the sonic point of 
the flow. 




1. INTRODUCTION 

Models of magnetohydrodynamic (MHD) disk accretion 
to a rotating star with an aligned dipole magnetic field 
have been discussed by many authors (e.g., Ghosh & Lamb 
1979; Camenzind 1990; Konigl 1991; Lovelace, Romanova, 
& Bisnovatyi-Kogan 1995). Major uncertainties remain 
however because of the many assumptions required to solve 
the MHD equations. For example, the validity of the 
standard formula for the disk's Alfven radius (Shapiro & 
Teukolsky 1983) is unknown and in fact this radius may de- 
pend on the rotation rate of the star (Lovelace, Romanova, 
& Bisnovatyi-Kogan 1999). 

Definite predictions can be made about funnel flows 
from disks to magnetized stars assuming a dipole mag- 
netic field. These flows shown in Figure 1 naturally ex- 
plain the large infall velocities observed in many T Tauri 
stars (Hartmann, Hewett, & Calvet 1994; Martin 1996; 
Najita et al. 2000). Such flows may also occur in disk 
accretion to rotating magnetized neutron stars (Li & Wil- 
son 1999). Here, a study is made of stationary, axisym- 
metric ideal magnetohydrodynamic (MHD) flows from an 
accretion disk to a rotating star with an aligned dipole 
magnetic field. The MHD equations lead to a Bernoulli 
equation which we analyze in the general case where pres- 
sure and density of a fluid particle are related as p cx p 1 
at constant entropy. A Bernoulli equation for the isother- 
mal case p tx p was discussed earlier by Li and Wilson 
(1999). Detailed calculations of emission line shapes were 
made by Hartmann et al. (1994) and Martin (1996) for 
a model where plasma free-falls along dipole field lines of 
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Fig. 1. — The thefigure shows a plasma accretion disk corotat- 
ing with the aligned magnetic dipole field of a star. Plasma flows 
from the disk to the star along a flux tube of the dipole, going 
from subsonic to supersonic flow. Conservation of magnetic flux 
gives BpAp = const, and conservation of mass gives pv p A p = const, 
where v p is the plasma speed along the field line and A p is the 
cross sectional area of the flux tube. The dashed line marked with 
(*) indicates the sonic point and that marked with (s) indicates a 
standing shock. 

We use inertial, cylindrical (r, (f>, z) and spherical 
(R, 6, (j>) coordinate systems. As in the case of matter 
outflow along open magnetic field lines threading a disk 
(Ustyugova et al. 1999), there are a number of integrals 
of the plasma motion which are conserved along mag- 
netic field lines labeled by the flux function &(r, z) — 
rAcj, =const, where is the vector potential. The con- 
stant K (\&) arises from the conservation of mass, f2(^) 
from conservation of hclicity, A(\&) from conservation of 
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angular momentum, S(^) from conservation of specific en- 
tropy, and E($>) from conservation of energy (Bernoulli's 
constant) (see for example Ustyugova et al. 1999). 

v p = ^b p , n ( *) = w -^, (i) 

A(*)= W r 2 -^, S = S(*), (2) 



if 



S(*) = | + i(o;. 



2„2 



47rpr 
fl 2 r 2 



(3) 



where the p subscript indicates the poloidal (r, z) [or 
(i?, 0)] part of the vector, oj = v^/r the angular velocity of 
the matter, w is the specific enthalpy, and <f> g = —GM/R 
the gravitational potential of the star. (The mass of the 
disk is assumed negligible.) An additional equation - the 
Grad-Shafranov equation - follows from the equilibrium of 
forces across magnetic field lines (Lovelace et al. 1986). 

In §2 we discuss simplifications which occur in the strong 
field limit, and in §3 we analyze the Bernoulli equation (3). 
In §4 dimensionless variables are introduced. In §5 we dis- 
cuss the nature of the transonic solutions which go from 
subsonic flow near the disk to supersonic flow near the 
star and show that the behavior of the flows is different 
for 7 < 7/5 and 7 > 7/5. In §6, we describe the nature of 
the possible subsonic outflows from the disk to the star. In 
§7, we contrast the considered flows with spherical Bondi 
accretion. In §8 we discuss conditions for the strong field 
approximation to hold. In §9 we apply the theory to ac- 
cretion to a neutron star and to a young stellar object. In 
§10 we summarize the work. 

2. STRONG FIELD APPROXIMATION 

We investigate the case where the magnetic field of the 
star is dominant. Namely, we suppose that the star has 
dipolc magnetic moment n in the direction coinciding 
with rotation axis of the disk. Thus the poloidal mag- 
netic field in the considered region is the dipole magnetic, 
Bp = 3R(/x- R)/i? 5 - fi/R 3 . Further, we suppose that the 
magnetic field energy-density is dominant compared with 
matter energy-density, 

^ » ( p , pv 2 , p$ ) . (4) 

Hence the matter flow does not disturb the dipole mag- 
netic field significantly, and the toroidal component of the 
magnetic field is small, \B$\ <C |B P |. Consequently, the 
full Grad-Shafranov equation is not needed. In §8 we give 
necessary conditions for the strong field approximation to 
hold. 

In spherical coordinates, the flux function of the dipole 
magnetic field is ^ = p sin 2 (6*)/i?. Thus, the equation for 
the dipole field lines is 

R = # d (*)sin 2 0, (5) 
where Rdi't') = /V* i s the radius at which the mag- 
netic field line <J/ =const passes through the disk (z = 0). 
Note that B R = (R 2 sin6)- 1 d^/de = 2pcos9/R 3 , and 
B e = -(Rsmey^^d^/dR = pamO/R 3 . Figure 1 shows 
the envisioned geometry. 

Equations (1) and (2) can be solved for uo and B$ to 
give 

l-( PA /p)(h/r 2 ) _ n l^ h /r 2 



LU = il 



Pa/p v 1 



Pa/p ' 
(6) 



where pa = K 2 /{Air) and h = A/O. Note that pa/ P = 
V p/ V A P 1S ^e poloidal Alven Mach number squared calcu- 
lated using the Alfven velocity va p = \H p \/\/4np. Owing 
to our assumptions, pa/ P *C 1 everywhere in the consid- 
ered region. Consequently, 

(7) 

We assume further that {p A /p)\l — h/r 2 \ <C 1 with the 
result that \u> — Cl\ <C ft- From equation (13) note that 
|B*|/|B P | = (rn/v Ap )\l - h/r 2 \ « 1. 

Thus a given fluid particle moves along the poloidal 
magnetic field v p cx B p , and its angular velocity fl^) 
is a constant. The value S7(\I/) is the angular velocity 
of rotation of the footpoint of the magnetic field line 
"J/ =const, which we suppose is frozen into the star. Thus, 
Sl('I') = flgtar =const is the angular rotation rate of the 
star. We omit the subscript on in the following. 

3. Bernoulli's constant 

Consider now the matter flow along a specific magnetic 
field line R = Rd(^) sin 2 6. The effective potential along 
this line in a reference frame rotating with rate f2 is 

* ^ GM tt 2 R 2 sm 2 6 GM fl 2 R 3 



R 



R 2R d 



(8) 

The magnitude of the magnetic field along this field line is 
iUR) = ^(4-3sin 2 AV' 2 " 



3RY 



P_ 

R 3 ^ R? 
Thus the Bernoulli constant along this field line is 

E = F(R,p) ( - 3R 



— ) +w + <P(R) 



(9) 



(10) 
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where w — Sp 1 ^ 1 / (7 — 1) = c 2 /(j — 1), with c s the sound 
speed and S = "fp/ p 1 the specific entropy. 

We consider flows which are subsonic near the disk, 
I Up I *C c s for R < Rd- There are then three possible 
cases: (1) The flow is transonic, going from subsonic to 
supersonic at some distances i?* < R d ; (2) The flow re- 
mains subsonic between the disk and the star; or (3) no 
stationary solution to the Bernoulli equation exists. 

For the transonic flows, note that the p derivative of 
equation (10) is 

f =^£, (id 

d P R P 

in that (dw/dp) s = c 2 / p. In general F(R,p) has a mini- 
mum as a function of p at, say, p m (R), and E = F(R,p) 
has two solutions for p. The larger p solution has dFj dp > 
and corresponds to subsonic flow. The smaller p solution 
has dF/dp < and corresponds to supersonic flow. 

For the transonic flows, the conditions for a smooth 
transition through the sonic point are 



dF_ 

dR 



= 



dF 
~dp~ 



= 



(12) 



In order to pass from subsonic to supersonic flow, this 
extremum of F(R, p) must be a saddle point which cor- 
responds to Frr*F pp * — F 2 Rtf < 0, where the subscripts 
indicate partial derivatives. The conditions for this are 
given in §5. 
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4. DIMENSIONLESS VARIABLES 



With R measured in units of Rd, we can write the effec- 
tive potential as 



(R/Rd 



\R/R d 2a 2 
where a = r c /Rd, and the characteristic length 



(GM\ 



1/3 



/ M 
1.5 x 10 8 cm 



1/2 



p \ 2/3 

3 (13 > 



is the "centrifugal radius" where the centrifugal force bal- 
ances the gravitational attraction of the star ( in the z = 
plane in the absence of magnetic forces), and P — 27r/f2 is 
the star's rotation period. For a young star of solar mass, 
r c w 2.9 x 10 n cm(P/ld) 2 / 3 . 

Useful dimcnsionlcss variables are 



R = 



E = 



R_ 

R~ d ' 

E 



„ 4iranR* 
P= — T? P 



5 = 



Kp, 

s (- 



Kp 



fl 2 r 2 \ATtaVLR d 



(14) 



and F = F/(nr c ) 2 , (c s ) 2 



Sp 



v p /(nr c ) = (4 - 3R) 1/2 /(pR 3 ), where a = r c /R d . Note 
that at R = 1, /w p = 1. For simplicity of notation we now 
drop the hats. We then have R s tar/Rd < R < 1, where 
Rstar is the star's radius. 

In terms of these variables, Bernoulli's equation becomes 



E = F(R,p) = 



4-3i? Sp 
~2p 2 W 1 



+ 



7-1 



R + 2a 2 



(15) 



We suppose that the disk matter is at a relatively low tem- 
perature, that is, w(R d ) = c 2 s (R d )/(j- 1) < -$(i? d ), and 
that the outflow speed from the disk is small v 2 (Rd) <C 
— <&(R d ). Under these conditions the dimensionless quan- 
tities at the starting point of the flow (R = 1) are: 
the Bernoulli constant, E « —a — l/2a 2 , the sound 
speed, c s (l) <C 1, and the density, p(l) » 1. Because 
5 = c 2 /p~<- 1 , we also have 5 = 5(1) < 1. 

5. TRANSONIC FLOWS 

The mentioned conditions for a smooth transition 
through the sonic point give 



= 







of 

dR 
OF 
~dp~ 



3(8 -5i?*) , a 
2plRj 

pIR* 6 



+ s P r 2 



3i?* 
2a 2 



(16) 
(17) 



In general, F pp * > 0. In order for this critical point to be a 
saddle point of F(R, p), we must have Frr* < F 2 Rtf /F pp ^. 
This requirement is satisfied if R* > (N/ D) 1 / i {2a 3 /3) 1 / 4 , 
where TV = WjR 2 - 45i? 2 + 140.R* - IOO7.R* + 8O7 - 112 
and D = m^Rl - 15i? 2 + 52i?, - 188 7 i?* + 144 7 - 48. This 
condition on R* is always satisfied if 7 < 7/4. For 7 = 5/3, 
the condition is satisfied for R* > 0.581(2a 3 /3) 1 / 4 . The 
nature of the function F(R, p) can be understood from its 
contour plot shown in Figure 2a for a sample case. 




Fig. 2. — Contours of the Bernoulli function E = F(R,p) in 
the top panel (a) and E = F(R,v p ) in the bottom panel (b), for 
a sample case with 7 = 5/3, dimcnsionlcss entropy S = 0.005, and 
o = l. Distance from the center of the star R is normalized by R^. 
The transonic flow starts at R = 1 and follows the dashed line. The 
sonic point where v p = c s - indicated by the asterisk - is the saddle 
point of the surface F at (R*,p*). For R > the flow is subsonic 
and dF/dp > 0, while for the reversed inequalities the flow is super- 
sonic. A possible subsonic flow would follow the upper E = —1.42 
contour in (a) and the lower E = — 1.42 contour in (b). 

The term in equation (16) involving can be rewrit- 
ten in terms of 5 using equation (17). This term is then 
seen to be negligible compared with the terms involving a 
because 5 <C 1. Equation (16) is then seen to correspond 
to a balance of the gravitational attraction and centrifugal 
force along the field line. Under this condition, equation 
(16) can be satisfied only for 

a < a = (3/2) 1/3 w 1.145 . 

Consequently, 



(18) 



R* ~ — tt~ 



2a 3 
3 



1/4 



s~— 



4:-3R ^\-<+ 1 



4 - 3iL \ 



(19) 



Thus the mass flux-density is = /0*c s * = (4 — 

3-R,) 1 / 2 / R 3 . The approximation involves neglecting the 
p„ term in equation (16) which is valid for 5 <C 1 and a 
not too small compared with unity. In the next approx- 
imation, we find R* « (2a 3 /3) 1 / 4 [l - 1.15 3 / 4 a - n / 8 ] for 
7 = 5/3, and R* « (2a 3 /3) 1 / 4 !! - 0.885°- 87 a -°- 84 ] for 
7 = 1.3, and both expressions for i?» <C 1. If the fac- 
tor in the square brackets is less than about 0.863 then 
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as mentioned below equation (17) the saddle point disap- 
pears (by becoming a minimum). The saddle point given 
by equation (19) is the only saddle point with R* < 1. 
Figure 3 shows the dependences of i?* and c s * obtained 
by solving equations (16) and (17) numerically without 
approximations . 
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Fig. 3. — Dependences of the sonic point distance i?* and sound 
speed c* on a = r c /R d for the cases of 7 = 5/3 and 7 = 1.3 and 
dimensionless entropy S = 0.005. The dashed line is for the approx- 
imation of equation (19). For the 7 = 5/3 case, the saddle point of 
F(R, p) changes into a minimum for a < 0.218 so that these values 
are not of interest. 

Figures (4) and (5) show sample radial profiles of the 
fluid variables for 7 = 5/3 and 1.3, respectively. A 
value of 7 smaller than the ideal monotonic gas value 
is of interest for a high temperature plasmas where the 
electron heat conduction can be modeled by smaller 7 
(as is the case for the solar wind). The difference 
in behavior of the cases shown in Figures (4) and (5) 
can be understood by considering the Bernoulli equa- 
tion for R < 1. In this limit v p w (2a/i?) 1/2 [in 
dimensional variables this is v v « [2GM/R) 1 / 2 , which 
is the free-fall speed], p « (2/a) 1 / 2 R^ 5 / 2 , and c s f» 
{2/a)^~ 1 ^ i S 1 l 2 R~ b ^' 1 ^ i . Hence the Mach number 
varies as M = v p /c s « 2(a/2)(''+ 1 ^ 4 S~ 1 / 2 R^- r ^ 1 . 
Thus for 7 > 7/5, the Mach number initially increases with 
decreasing R, but for small enough R, M. decreases with 
decreasing R. The radius of the maximum of M. is found 
tobeR max m (5 7 -7)/{6(7-l)[l + l/(2a 3 )]}. The maxi- 
mum value of M. is therefore oc S^ 1 / 2 . On the other hand 
for 7 < 7/5, A4 increases monotonically with decreasing 
R. Therefore, the changeover in behavior between Figures 
(4) and (5) occurs at 7 = 7/5. The radial variation of 
both v p and M. are clearly important for calculations of 
the emission line profiles (Hartmann et al. 1994). 

For 7 > 7/5, the Mach number does not decrease 
through M. = 1. For R <C 1, Bernoulli's equation is 
2/(p 2 R e ) + S'p 7 ~ 1 /(7 - 1) « a/R. Substituting p = r]/RP 
with p = 6/(7 + 1) gives 2/rj 2 + Sy~7(7 ~ 1 ) = L (v) ~ 
aJ?( 5 T _7 )/(i' _1 ). Clearly, L(rf) has a minimum value of 
L min = [2 + 4/(7- l)](S/4) 2/(7+1) at Vmm = (4/5) 1 /(7+D 
which corresponds to M. = 1. Thus for 7 > 7/5, there is 
a minimum radius where a stationary flow exists, R m in = 
[(2/a)(7 + l)/(7 - l)](T+ 1 )/( 5 'r- 7 )(5/4) 2 /( 57 - 7 ). For most 
conditions this radius is expected to be less than the radius 
of the star. 

Near the star the flow is expected to have a standing 
shock where M. discontinuously decreases to a value less 
than unity. For 7 < 7/5 there also must be a shock. Note 
that the ratio of the temperatures across this shock for an 
ideal gas is T 2 jT x = 1 + 2(7- 1)( 7 7W 2 + 1){M 2 - l)/[(7 + 
l) 2 M 2 ] (w 27(7 - l)A4 2 /( 7 + l) 2 for M > 1), where T x 
is the upstream temperature and T2 is downstream (closer 



to the star). Calculation of the location of this shock is 
beyond the scope of this work. 
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Fig. 4. — The thefigure shows the radial variations of the 
dimensionless density p, poloidal velocity v p , and Mach number 
M = Vp/c s for cases with 7 = 5/3, a = 1, and dimensionless 
entropies S = 0.001 (top panel, a) and S = 0.005 (bottom panel, 
b). For (a), E fa -1.462, c si fa 0.150, p d 129.4, R, ss 0.898, 
and c s , fa 0.0841. For (b), E fa -1.428, c sd fa 0.218, p d fa 29.4, 
R-t fa 0.884, and c s * fa 0.156. The dashed curves represent 
constij - 5//2 fitted to a point at the top of the plot. 
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Fig. 5. — Radial variations of the dimensionless density p, poloidal 
velocity v p , and Mach number M = v p /c s for a case with 7 = 1.3, 
o = l, and dimensionless entropy 5 = 0.005. The dashed curve 
represents constiJ _5//2 fitted to a point at the top of the plot. 
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Fig. 6. — Dimensionless function in equation (20). The dimen- 
sionless sound speed c s( j must be less than C(a). 

The value of entropy S for transonic flows St s is deter- 
mined by conditions at origin of the flow at the disk where 
v pd ^ C L> w here c sd is the sound speed at the disk. At 
the disk. 



E 



7-1 



1 

2c^ 



while at the sonic point, 



17 + 1 2 _{a_,]$_ 

2 7 -l C " \R, + 2a 2 

Note that c 2 „ = Spl' 1 and p, = (4 - 3JZ.) 1/2 /(c«-H2)- 
Therefore, in order to have a transonic flow, we need 



S — St.* = 



7 + 1/ \4-3R 



ir 



[c 2 sd ~C 2 (a) 



where 



(20) 



C 2 (a) = ( 7 -1) [a+^-aa 
la 



is obtained by using equation (19), and c\ — (3/2) 1 / 4 + 
(2/3) 3 / 4 /2 w 1.476. Figure 6 shows C(a). Clearly, for a 
given value of a it is necessary to have sufficiently high 
disk temperature or disk sound speed, c sd > C(a). This 
inequality and inequality (18) for transonic flow can be 
satisfied only for 



< a < ao 



(21) 



where a m i n is such that c s d = C{a m i n ). Inequalities 
(21) correspond to a limited range of Rd, r c /aa < Rd < 
f c /a m in- For example, for constant c s d = 0.1 and 7 = 5/3, 
we have c s d > C(a) for 1.03 < a < 1.145 and S varies be- 
tween and 0.00147; for c s d — 0.2, the allowed range is 
0.928 < a < 1.145 and S varies between and 0.00933. 



6. SUBSONIC FLOWS 

For a m i n < a < ao (with a m i n given above), it is 
easy to show that there are outflows from the disk to the 
star which are subsonic. Figure 7a shows the R depen- 
dences for such a case. For E« 1, equation (15) gives 
c s « [a(7 - 1)/ R] 1 ' 2 and p « [0(7 - ^/(SR)] 1 '^-^. 

For a < a m i n , there are no stationary solutions to the 
Bernoulli equation (15). Equivalently, there are no sta- 
tionary outflows for Rd > r c /a m i n . 

For a > ao or Rd < r c /ao, there are also subsonic out- 
flows from the disk to the star as shown in Figure 7b. 
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Fig. 7. — Radial dependences of density p, sound speed c s , and 
Mach number for subsonic flows from disk to star. The top panel 

(a) has a = 1.0, S = 0.001, 7 = 5/3, and c s[ j fa 0.173. In this case 
the Mach number rises to about 0.26 at R = 0.9. The bottom panel 

(b) has a = 1.3, S = 0.001, 7 = 5/3, and c sd fa 0.1. 

7. CONTRAST WITH BONDI ACCRETION 

The Bernoulli equation(lO) differs from that for spher- 
ical Bondi (1952) accretion in that here the plasma flows 
along the magnetic field lines, the matter is rotating, and 
the flow starts from a finite radius of a disk. For the Bondi 
flow, mass conservation gives AirR 2 pvR = M =const so 
that the Bernoulli equation is 



E 



7-1 



_ (M/47T) 2 Sp 

2p 2 R* 7 - 1 
/pjT 1 , with c.c 



GM 
~R~ 



(22) 



where S = cJoo/pjT > with c soo and the sound speed 
and density at a large distance from the star. In con- 
trast, the first term of equation (10) is oc l/(p 2 R 6 ) ow- 
ing to the "focusing" of the magnetic field, and further, 
equation (10) includes the centrifugal potential — fi 2 r 2 /2. 
For small R, the Bondi flow has v R » (2GM/R) 1 / 2 as in 
the for the up and over flow; however, p ~ \jR?l 2 and 
c s ~ 1/i? 3 ^ 7-1 )/ 4 , whereas for the funnel flow p ~ 1/i? 5 / 2 
and c s - 1/i? 5 ^ 1 )/ 4 . Thus, the Mach number for the 
Bondi flow M - 1/R (5 - 3 ^/ 4 always increases (or is con- 
stant) with decreasing R for 7 < 5/3. 

8. CONDITIONS FOR STRONG FIELD APPROXIMATION 

In the strong field approximation of §2, the disk plasma 
is assumed to corotate with the star even though r is 
less than the corotation radius r c . Equilibrium of this 
part of the disk thus requires an additional outward ra- 
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dial force which arises naturally from the slight bend- 
ing of the field lines passing through the disk. This ra- 
dial magnetic force (per unit area of the disk) is simply 
prnag = -(B r B z ) h /2ir > 0, where the /i-subscript in- 
dicates evaluation at the disk surface z = h [equation 
(5) of Lovelace et al. 1995]. Because (r c — r)/r c <C 1, 
the required radial magnetic field is small, (B r /B z )h = 
(6iran 2 /B 2 z )(r c -r) = i{VLr c /v A ) 2 h{r c -r)/r 2 c < 1, where 
va is the Alfven velocity based on the midplane density of 
the disk and B z (r, 0) . 

Also, in the approximation of §2 the toroidal magnetic 
field is neglected. A small deviation of the disk rotation 
rate u> d (r) from the star's rate O can lead to a signifi- 
cant toroidal magnetic at the disk, (B^)h = —{hr /rf){ujd — 
ri)B Zl where r\ is the magnetic diffusivity of the disk [equa- 
tion (2) of Lovelace et al. (1995)]. If r\ is of the order of 
the viscosity v given by the Shakura and Sunyaev (1973) 
prescription v — ac s h (with a < 1), then, in order to have 
\{B$)h/B z \ <C 1, we must have \r(uj d - 0)| <C ac s . 

9. APPLICATIONS 

Table 1 summarizes the relevant parameters for the 
cases of disk accretion to (1) a rotating magnetized neu- 
tron star, and to (2) a rotating young stellar object. In 
both cases, the Alfven radius ta (Shapiro & Teukolsky 
1983) is larger than r c so that the disk plasma is expected 
to corotate with the star. The disk plasma may move 
radially across the magnetic field by the interchange in- 
stability (Kaisig, Tajima, & Lovelace 1991; Rastatter & 
Schindler 1999). For the neutron star the surface mag- 
netic field (fj,/R^ tar ) is assumed to be 10 12 G, while for the 
young star it is 10 4 G. 

For the neutron star case, we assume c s d/{Slr c ) = 0.1 
(with dimensions restored), which corresponds to a surface 
temperature of the disk Td ~ 4.1 x 10 7 K for a fully ionized 
hydrogen plasma and 7 = 1.3 with mean particle mass 
m#/2. The half-thickness of the disk is h ~ c s /Q, w 0.1r c . 
From Figure 6 with 7 = 1.3, we find that transonic out- 
flows from the disk to the star are possible for a in the 
range a — r c /Rd ~ 0.95 — 1.145. which corresponds 
to Rd ~ (0.873 — 1.053)r c . The sonic point location 
R*/Rd = sin 2 (6>*) has 0» varying between 69° (the lower 
limit of a and the upper limit of Rd) and 90° where the 
sonic point is inside the disk. The height of the sonic 
point z* = Rd sm 2 (9 if ) cos(0*) is < 0.273r c . Thus, tran- 
sonic outflows occur over an interval ARd ~ 0.18r c . In 
dimensionless variables, pdv p d = 1- Therefore, if all of 
the disk accretion goes into transonic flow along the field 
lines, M acc = 2(2ir)RdAR d K jt/ '(AirR^), where one of the 
factors of two in the numerator accounts for the two sides 
of the disk. Thus the characteristic density normalizing p 
in equation (14) is 

Kp _ M acc 

P ° ~ AiranRj _ AnaVLR 2 d AR d ' ( ' 

For the mentioned values and those in Table 1, we find 
po w 1.4 x 10 _9 g/cm 3 . The ratio of the maximum to 
minimum radii of the flow is r c /R star — 150 and that 
the free-fall speed at the star is v/f = (GM/Rstar) 1 ^ 2 ~ 
1.1 x 10 10 cm/s. The angle at which the flow hits the star 
is 6 st ar = sin _1 [(i? star /r c ) 1 / 2 ] w 4.7°. At this point the 
flow velocity is at an angle w 9 star /2 with respect to R. 



For the case of a young stellar object, we assume 
c s d/(Qr c ) — 0.05, which corresponds to a surface tem- 
perature of the disk T d w 2, 800K for a slightly ionized 
hydrogen plasma and 7 = 5/3. The half-thickness of the 
disk is h ~ c s /Q rj 0.05r c . The height of the sonic point is 

< 0.174. The considerations for this case are similar to 
those for the neutron star. The range of a is rs 1.08 — 1.145, 
R d « (0.873 - 0.926)r c and AR d ~ 0.053r c . The angle to 
the sonic point varies from 9* rj 78° (at the lower limit 
of a) to 90°. The characteristic density in this case is 
po ~ 10~ 12 g/cm 3 , assuming that all of the disk accretion 
goes into transonic flow along the field lines. The ratio 
r c / ' Rstar = 8.6 is much smaller than for the neutron star 
case, and Vff(R s t ar ) ~ 365 km/s. The angle at which the 
flow hits the star is 9 star ~ 20°. At this point the flow 
velocity is at an angle w 10.4° with respect to R. 

10. CONCLUSIONS 

This work analyzes the stationary ideal magnetohydro- 
dynamic funnel flows along the magnetic field lines of a 
rotating star with an aligned dipolc magnetic field. Isen- 
tropic flow is assumed so that a given plasma blob main- 
tains p/p 1 = const with 1 < 7 < 5/3, and the magnetic 
field is assumed strong, pv 2 <C B 2 /47r. Earlier, Li and 
Wilson (1999) considered the Bernoulli equation assuming 
p oc p. Close to the disk the flow is subsonic. There are 
three possible cases: (1) stationary transonic flow where a 
sonic point exists along the dipolc field line; (2) stationary 
subsonic flow; and (3) no stationary flow exists. Over a 
range of Rd < r c , subsonic inflows are possible, where R d 
is the radius at which the dipole field line intersects the 
disk and r c is the 'corotation radius.' For Rd >^ r c there 
are no stationary flows. 

The transonic flows, which become free-fall inflow near 
the star, are possible only for a narrow range of Rd ~ r c . 
At the sonic point these flows have an approximate bal- 
ance of the centrifugal and gravitational force along the 
field line. The flows have different behaviors for 7 > 7/5 
and < 7/5. For 7 > 7/5, the Mach number M = |v|/c s 
initially increases with R decreasing from R d , but for R 
decreasing from 0.22Rd (for 7 = 5/3) the Mach num- 
ber decreases. In the other case, 7 < 7/5, M. increases 
monotonically with decreasing R. The variation of both v p 
and M are clearly important to the nature of the standing 
shock near the star and to the determination of emission 
line profiles. We discuss numerical values for funnel flows 
to neutron stars and young stellar objects. For the neu- 
tron star case we argue that 7 is likely to be less than 
7/5. This combined with the large value of r c /r star im- 
plies that the Mach number of the flow approaching the 
star should be very much larger than unity. On the other 
hand for flows to a forming star we expect 7 > 7/5. This 
combined with the not very large values of r c /r sta r implies 
only moderately large Mach numbers near the star. The 
flow solutions analyzed in this work are also of interest for 
comparison with MHD simulations. 
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Table 1 
Sample Parameters' 1 



p 


Rstar 


M 


M 


Lace 


r c 


r A 






7 


1 s 
5 d 


10 6 cm 
10 n cm 


1M 
1M 


l(T 9 M Q /yr 
l(r 7 M /yr 


0.85xl0 37 erg/s 
0.85xl0 34 erg/s 


1.5xl0 8 cm 
8.6xlO n cm 


3.7xl0 8 cm 
9.6xlO n cm 


10 30 Gcm 3 
10 37 Gem 3 


9.4xl0 8 cm/s 
125 km/s 


1.3 

5/3 



a P = 27r/ri is the star's rotation period; R s tar its radius; M its mass; Al is the accretion rate; L — G M Af / v star 
is the accretion luminosity; r c = (GM/fi 2 ) 1 / 3 is the 'corotation radius'; ta = [p 4 /(2GMM 2 ))] 1 / 7 is the Alfven radius 
of the disk (see Shapiro and Teukolsky 1983); [i is the magnetic moment of the star; ilr c is the azimuthal velocity of 
the disk at the corotation radius; and 7 is the usual adiabatic exponent. 



